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Scalar and tensor glueball spectrum is studied using an improved gluonic action on asymmetric lattices in the 
pure SU(3) gauge theory. The smallest spatial lattice spacing is about 0.08/m which makes the extrapolation 
to the continuum limit more reliable. In particular, attention is paid to the scalar glueball mass which is known 
to have problems in the extrapolation. Converting our lattice results to physical units using the scale set by the 
static quark potential, we obtain the following results for the glueball masses: Mg(0 ++ ) = 1730(90)MeU for the 
scalar glueball mass and M G (2 ++ ) = 2400(95) AfeV for the tensor glueball. 



1. INTRODUCTION 

It is believed that QCD is the theory which 
describes strong interactions among quarks and 
gluons. A direct consequence of this is the ex- 
istence of excitations of pure gluonic degrees of 
freedom, i.e. glueballs. However, due to their 
non-perturbative nature, the spectrum of glue- 
balls can only be investigated reliably with non- 
purterbative methods like lattice QCD (l]-|6) . Re- 
cently, it has become clear that such a calculation 
can be performed on a relatively coarse lattice 
using an improved gluonic action on asymmet- 
ric lattices [Q-^]. In this paper, we present our 
results on lowest scalar and tensor glueball spec- 
trum calculation. The spatial lattice spacing in 
our simulations ranges from 0.08/m to 0.25/m 
which enables us to extrapolate more reliably to 
the continuum limit. The improved gluonic ac- 
tion we used is the tadpole improved gluonic ac- 
tion on asymmetric lattices as described in 
It is given by: 



In the above expression, (3 is related to the bare 
gauge coupling, £ = a s /a t is the (bare) aspect 
ratio of the asymmetric lattice with a s and o* 
being the lattice spacing in spatial and tempo- 
ral direction respectively. The parameter u s is 
the tadpole improvement parameter to be deter- 
mined self-consistently from the spatial plaque- 
ttes in the simulation. Py and Pqi are the spatial 
and temporal plaquette variables. i?y designates 
the 2x1 Wilson loop (2 in direction i and 1 in di- 
rection j). Using spatially coarse and temporally 
fine lattices helps to enhance signals in the glue- 
ball correlation functions. Therefore, the bare as- 
pect ratio is taken to be some value larger than 
one. In our simulation, we have used £ = 3 for our 
glueball calculation. It turns out that, using the 
non-perturbatively determined tadpole improve- 
ment 1 10 1 parameter u s , the renormalization ef- 
fects of the aspect ratio is small , typically of 



the order of a few percent for practical values of (3 
in the simulation, which is ignored in this paper. 
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2. MONTE CARLO SIMULATIONS 

We have utilized a Hybrid Monte Carlo algo- 
rithm to update gauge field configurations. Sev- 
eral lattice sizes have been simulated and the de- 
tailed information can be found in Table. 1. For 
each lattice with fixed bare parameters, order of 
a few thousand configurations have been accu- 
mulated. Each gauge field configuration is sep- 
arated from the previous one by several Hybrid 
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Table 1 

Simulation parameters and the corresponding lat- 
tice spacing in physical units obtained from Wil- 
son loop measurements. Parameters used in the 
smearing process for Wilson loop measurements 
are also listed. 



Lattices 





\w 




% 




ro/a s 


8 3 x 24 


2.4 


0.20- 


- 0.40 


2 ~ 


4 


1.98(2) 


8 3 x 24 


2.6 


0.20- 


- 0.40 


2 ~ 


4 


2.48(2) 


8 3 x 24 


3.0 


0.20- 


- 0.40 


4~ 


G 


4.11(4) 


8 3 x 24 


3.2 


0.25 - 


- 0.50 


4~ 


G 


5.89(8) 


10 3 x 30 


3.2 


0.25 - 


- 0.50 


4~ 


G 


5.89(8) 



Monte Carlo trajectories, typically 5 ~ 10, to 
make sure that they are sufficiently de-correlated. 
Further binning of the data has been performed 
and no noticeable remaining autocorrelation has 
been observed. 

2.1. Setting the scale 

In our simulation, the scale is set by measur- 
ing Wilson loops from which the static quark 
anti-quark potential V(R) is obtained. Using the 
static potential between quarks, we are able to 
determine the lattice spacing in physical units 
by measuring ro , a pure gluonic scale determined 
from the static potential pl]| . The definition of 
the scale r is given by: R 2 dV{R)/ dR\n =ra = 
1.65. In physical units, ro is roughly 0.5 fm which 
is determined by comparison with potential mod- 
els. For a recent determination of ro, please con- 
sult Ref. @. 

In order to measure the Wilson loops accu- 
rately, it is the standard procedure to perform 
single link smearing on the spatial links and 
to utilize thermally averaged temporal links of the 
configurations. The smearing scheme can be per- 
formed iteratively on the spatial links of gauge 
fields for as many as nw times with a smear- 
ing parameter \\y . The smearing parameters Aw 
and nw used in this process are also listed in Ta- 
ble. 1. Wilson loops are then constructed using 
these smeared spatial links and thermally aver- 
aged temporal links. For a Wilson loop of size 



R x T, it is fitted against: 

W(R,T) T ~°° Z(R)e- v W T . (2) 

The static quark potential V(R) is obtained from 
the effective mass plateau in temporal direction. 
Wilson loops measured along different lattice axis 
and at different lattice spacings are seen to lie on 
a universal line which is an indication that the 
improved action restores the rotational symmetry 
quite well. 

The static quark potential is fitted according 
to a Coulomb term plus a linear confining poten- 
tial which is known to work well at these lattice 
spacings |9| . The potential is parameterized as: 

V(R) = V + e c /R + aR . (3) 

From this and the definition of ro, it follows that: 

/l.65 + e c 

To convert the measured result to ro /a s , we have 
also used the value of £ taken as the bare value. 
The results of the spatial lattice spacing in physi- 
cal units are also included in Table. 1. The errors 
for the ratio ro/a s are obtained by blocking the 
whole data set into smaller blocks and extracting 
the error from different blocks. 

2.2. Glueball mass measurement 

To obtain glueball mass values, it is necessary 
to constructed glueball operators in various sym- 
metry sectors of interests. Scalar glueball is in 
the Ai representation of the cubic group; tensor 
glueball is in representation E + T2, which forms 
a 5-dimcnsional representation; vector glueballs 
are in the representation T\ |ll|,|3| . 

Glueball correlation functions are noisy and 
difficult to measure. In order to enhance the 
signal of glueball correlation functions, smear- 
ing and fuzzying have to be performed on spa- 
tial links of the gauge fields @,H,|). These tech- 
niques greatly enhance the overlap of the glue- 
ball states and thus provide possibility of mea- 
suring the mass values. The smearing and fuzzy- 
ing process involves typically 6—12 single link 
smearings plus a double link fuzzying. In our 
simulations, after performing single link smearing 
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Figure 1. The Wilson loop shapes used in con- 
structing the glueball operators. 



and double link fuzzying on spatial links, we first 
construct raw operators at a given time slice t: 
= Ex w x*t}> where W^l are closed Wil- 
son loops originates from a given lattice point 
x = (t,x). The loop-shapes studied in this calcu- 
lation are shown in Fig.l. Then, elements from 
the cubic group is applied to these raw operators 
and the resulting set of loops now forms a ba- 
sis for a representation of the cubic group. Suit- 
able linear combinations of these operators are 
constructed to form a basis for a particular ir- 
reducible representation of interest ||. We de- 
note these operators as {Oa (*)} where R labels 
a specific irreducible representation and a labels 
different operators at a given time slice t. 

In order to maximize the overlap with one 
glueball state, we construct a glueball opera- 



v, 



( a n >(t), where 6 {R) (t) = 
The coefficients 



tor 0(«)(t) = £ c 
O (R) (t)-(0\O (R) (t)\0 
are determined from a variational calculation. To 
do this, we construct the correlation matrix: 



C o/s (t) = X;<0|OW(t + r)OW(r)|0 



(5) 



The coefficients Da are chosen such that they 
minimize the effective mass 

(«)„,(#), 



(tc) 



1 



log 



V 







>C a0 (t c ) 



(R) /„, 
V" a 'Vp C a /3(0) 



(6) 



where tc is time separation for the optimization 
and repeated indices are summed over. In our 



simulation tc = 1 is taken. If we denote the op- 
timal values of by a column vector , this 
minimization is equivalent to the following eigen- 
value problem: 



r {R) 



(7) 



1R) 



with the lowest effective 

(R) 



The eigenvector Vq 

mass then yields the coefficients Vq™' for the oper- 
ator Qq R ^ (t) which best overlaps the lowest lying 
glueball in the channel with symmetry R. Higher- 
mass eigenvectors of this equation will then over- 
lap predominantly with excited glueball states of 
a given symmetry channel. 

With these techniques, the glueball mass values 
are obtained in lattice units and the final results 



Table 2 

Glueball mass estimates for the symmetry chan- 
nel A± + , E ++ and at various lattice spac- 
ings. The entries corresponding to the highest 
j3 value are the values after the infinite volume 
extrapolation. The last row tabulated the con- 
tinuum extrapolated result of the glueball mass 
values in units of 1/tq. 



(3 a t M A++ a t M_ 



E++ 



a t M, 



2.4 0.552(8) 

2.6 0.482(12) 

3.0 0.322(8) 

3.2 0.233(7) 



A 2- 



0.980(10) 1.002(9) 

0.760(16) 0.798(15) 

0.460(13) 0.470(13) 

0.323(12) 0.340(10) 



oo 4.23(22) 5.77(34) 5.92(32) 



are listed in Table. 2. The errors are obtained by 
binning the total data sets into several blocks and 
doing jackknifc on them. 

2.3. Extrapolation to the continuum limit 

Finite volume errors are eliminated by perform- 
ing simulations at the same lattice spacing but 
different physical volumes. This also helps to 
purge away the possible toleron states whose en- 
ergy are sensitive to the size of the volume. A sim- 
ulation at a larger volume is done for the small- 
est lattice spacing in our calculation. We found 
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that the mass of the scalar glueball remains un- 
changed when the size of the volume is increased. 
The mass of the tensor glueball is affected, which 
is consistent with the known result that tensor 
glueballs have a rather large size and therefore 
feel the finitencss of the volume more heavily. The 
infinite volume is obtained by extrapolating the 
finite volume results using the relation Q : 

at M(«) = a t M^ (l - ^ exp(-^ , (8) 

(A + + ) 

where z = 1 L s . Using the results for the 
mass of the E++ and T 2 glueballs on 8 3 24 and 
10 3 30 lattices for the same value of (3, the final 
result for the mass of these glueball states are 
obtained. Glueball mass values for other symme- 
try sectors are not sensitive to the finite volume 
effects. Therefore, in Table. 2, only the extrapo- 
lated values for the smallest physical volume are 
tabulated. Other entries are obtained from 8 3 24 
lattice results. 

As for the finite lattice spacing errors, spe- 
cial attention is paid to the scalar glueball sec- 
tor where the continuum limit extrapolation was 
known to have problems. Due to the simulation 
points at small lattice spacings, around 0.1/m 
and below, the ambiguity in this extrapolation 
is reduced. We have tried to extrapolate the re- 
sult using different formula suggested in Ref. ||], 
the extrapolated results are all consistent within 
statistical errors. For definiteness, we take the 
simple form: 



M G (a s ) = M G (0) + 



cr 








ro 









,(9) 



and the result is illustrated in Fig. 2. The final ex- 
trapolated results for the glueball mass values are 
also listed in Table. 2. The data points from our 
simulation results are shown with solid symbols 
and the corresponding extrapolations are plotted 
as solid lines. It is also noticed that the extrap- 
olated mass values for E ++ and T" 2 ++ channels 
coincide within statistical errors, indicating that 
in the continuum limit, they form the tensor rep- 
resentation of the rotational group. For compari- 
son, corresponding results from || are also shown 
with open symbols and the dashed lines. If one 
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Figure 2. The continuum limit extrapolation of 
glueball mass values in scalar and tensor channels. 
The solid symbols are results from this calculation 
with the corresponding continuum limit extrapo- 
lation represented by the solid lines. For compar- 
ison, the corresponding results from j^J^] are also 
shown with open symbols and dashed lines. 



would extrapolate linearly in (a s /ro) 2 using three 
data points with smallest lattice spacing, the re- 
sults are statistically consistent with the results 
using the extrapolation (Q) within errors. Our 
final extrapolated result for the scalar glueball 
mass lies higher, though still statistically consis- 
tent, than that of Ref.7. This has to do with 
the fact that simulation results at small lattice 
spacings are higher than those at larger lattice 
spacings. If one would extrapolate linearly using 
only the two data points with the smallest lattice 
spacing from one would arrive at a higher 

scalar glueball mass value and closer to our re- 
sult. This indicates that the difference between 
our final result on the scalar glueball mass and 
that of basically comes from the ambiguity of 
the extrapolation when one tries to include more 
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data points from larger lattice spacings. In our 
calculation, with the help of one more point at 
smaller lattice spacing, using three data points 
with smallest lattice spacings, one can extrapo- 
late linearly towards the continuum and obtain a 
result consistent with the result using a quadratic 
extrapolation. In the tensor channel, a constant 
extrapolation will also yield a consistent result. 
It is seen that, due to data points at lattice spac- 
ings around 0.1 fm and below, the uncertainties 
in the extrapolation for the glueball mass values 
are reduced. 

Finally, to convert our simulation results on 
glueball masses into physical units, we use the re- 
sult r^ 1 = AlOMeV. The errors for the hadronic 
scale ro is neglected. For the scalar glueball we 
obtain M G (0++) = 1730(90)MeV. For the ten- 
sor glueball mass in the continuum, we combine 
the results for the T%~ + and E ++ channels and 
obtain M G (2++) = 2400(95)MeV for the tensor 
glueball mass. 

3. DISCUSSIONS AND CONCLUSIONS 

We have studied the glueball spectrum at zero 
momentum in the pure SU (3) gauge theory using 
Monte Carlo simulations on asymmetric lattices 
with the lattice spacing in the spatial directions 
ranging from 0.08/to to 0.25/m. This helps to 
make extrapolations to the continuum limit with 
more confidence for the scalar and tensor glueball 
states. The mass values of the glueballs are con- 
verted to physical units in terms of the hadronic 
scale ro- We obtain the mass for the scalar 
glueball and tensor glueball to be: Mg(0 ++ ) = 
1730(90)M e y andM G (2++) = 2400(95)MeV. It 
is interesting to note that, around these two mass 
values, experimental glueball candidates exist. Of 
course, in order to compare with the experiments 
other issues like the mixing effects and the effects 
of quenching have to be studied. 
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